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Gravitational leasing as a possible explanation for some 
unidentified gamma-ray sources at high latitudes 

Diego F. Torres^ Gustavo E. Romero^ and Ernesto F. Eiroa ^ 

ABSTRACT 

We propose that some of the high-latitude unidentihed EGRET 7 -ray sources 
could be the result of gravitational lensing amplihcation of the innermost regions 
of distant, faint, active galactic nuclei. These objects have 7 -ray emitting regions 
small enough as to be affected by microlensing of stars in interposed galaxies. We 
compute the gravitational amplihcation taking into account effects of the host 
galaxy of the lens and prove that, whereas the innermost 7 -ray regions can be 
magnihed up to thousand times, there is no amplihcation at radio frequencies, 
leading to the observed absence of strong counterparts. Some new ehects in 
the spectral evolution of a gravitational microlensed 7 -ray AGN are predicted. 
Within a reasonable range of lensing parameters, and/or types of sources, both 
variable and non-variable EGRET detections at high latitudes can be explained 
by microlensing. The same phenomenon could also have an important incidence 
among the future GLAST detections at high-latitudes. 

Subject headings: gamma-rays: observations - galaxies: active galactic nuclei - 
cosmology: gravitational lensing 


1. Introduction 

The Third EGRET Gatalog of 7 -ray sources includes observations carried out between 
April 22, 1991 and October 3, 1995, and lists 271 point-like detections (Hartman et al. 
1999). About two thirds of them have no conclusive counterparts at lower frequencies, and 
are dubbed unidentihed. 
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These unidentified 7 -ray sources can be divided in two broad groups. The first one, at 
low latitudes, is probably related to several galactic populations such as radio-quiet pulsars, 
interacting supernova remnants, early-type stars with strong stellar winds. X-ray binaries, 
etc (see Romero 2001 and references therein). The second group of unidentified sources is 
formed by mid- and high-latitude detections. Gehrels et ah (2000) have shown -in a model 
independent way- that the mid-latitude sources are different from the bright population of 
unidentified sources along the galactic plane. Some of the mid-latitude detections (5° < | 6 | < 
30°) are thought to be associated with the Gould Belt (Grenier 2000, Gehrels et al. 2000), 
a starburst region lying at ~ 600 pc from Earth. A few other sources, at higher latitudes, 
might be the result of electrons being accelerated at the shock waves of forming clusters of 
galaxies (Totani et al. 2001). However, for most of the high-latitude unidentified sources 
(HL-UnidS), no other explanation seems to be available than they are AGNs yet unnoticed 
at lower energies (Reimer & Thompson 2001). 

A 7 -ray AGN population consisting of 66 members has already been detected by EGRET 
(Hartman et al. 1999). Photon spectra and variability indices of the HL-UnidS match well 
those of known 7 -ray AGNs. In Figure 1, we show the corresponding distributions, with the 
AGN-histogram adapted from Torres et al. (2001a). The variability criterion used here is the 
/-scheme, but other variability criteria yield similar results (Torres et al. 2001b). Population 
studies have already remarked that part of the sample of HL-UnidS is consistent with an 
isotropic population, a fact that also supports an extragalactic origin (Ozel & Thompson 
1996). 

All identified 7 -ray AGNs are also strong radio sources with fiat spectrum, as expected 
from synchrotron jet-like sources where the 7 -ray flux is the result of inverse Gompton 
scattering (Mattox et al. 1997). We could ask, then, why the HL-UnidS are not detected at 
lower frequencies, particularly in the radio band, if they are also AGNs? Here we propose 
that some of these sources are the result of gravitational lensing amplification of background, 
high-redshift, active galactic nuclei; blazars whose 7 -ray emitting regions are small enough 
as to be affected by microlensing by stars in interposed galaxies. 


2. Gravitational lensing of active galactic nuclei 

Since AGNs have emission regions of different size for different wavelengths, we expect a 
differential magnification of light. The innermost regions of AGNs, responsible for the 7 -ray 
emission, have linear sizes x ~ lO^^-lO^® cm (e.g. Blandford and Levinson 1995). Significant 
gravitational lensing requires that Doi/Dos x x/2 < Re, where Doi and Dos are the angular- 
diameter distances between the lens ( 1 ), or the source (s), and the observer (o), and Re is 



the Einstein radins of the lens. For typical redshifts (say, ~ few tenths for the lens, and 
2 : ~ 0.8 — 3 for the sonrce), in a standard cosmological model, we hnd that stars of few 
solar masses have Einstein radii ~500 times bigger than the projected AGN’s 7 -ray emitting 
region, and thus, that the latter can be perfectly amplihed. At the same time, since radio 
emission is originated far down the jet, the sizes of the corresponding radio regions {x > 10 ^^ 
cm) exceed the Einstein radius of the lenses and lead to the absence of radio counterparts. 
In the case of the optical emission, if it is the same particle population giving rise -in the 
same region- to gamma-rays (through inverse Compton scattering) and optical emission 
(through synchrotron radiation), then we can also expect gravitational magnihcation of the 
optical luminosity. In some cases with high (lensing) optical depth, we could even expect 
simultaneous intensity variations in both bands, something that could in principle be tested 
with the improved capabilities and source location accuracy of the forthcoming GLAST 
satellite. Optical monitoring of stellar-like sources inside GLAST error boxes could lead 
to the identification of new gamma-ray blazars, and to some knowledge of their redshift 
distribution. One interesting fact is that the spectral evolution that we predict below will 
only happen in the gamma-ray band, since it critically depends on the size of the gamma-ray 
spheres (given by the opacity to pair production), a process not operative in the optical band. 
The idea that some extreme properties of distant AGNs can be the result of gravitational 
microlensing is not new (e.g. Ostriker & Vietri 1985, Gopal-Krishna & Subramanian 1991, 
Romero et ah 1995). However, microlensing effects upon the 7 -ray emission have not been 
discussed in detail yet. 

It is usually assumed that the a priori probability of Ending a small group of distant, 
gravitationally magnified objects is below 1%. Indeed, recent results (E.L. Turner, personal 
communication, 2001 ) taking into account the clustering of lenses in interposed galaxies, give 
values between 10 “^ - 10 “^ for the a priori probability of finding gravitational magnified 
sources in random directions of the sky. In addition, for those directions where there is 
gravitational lensing, the probability of it happening with optical depth above 0.2 is very 
high. In our case, the number of potential compact 7 -ray emitting background sources is 
extremely large: only the last version of the Veron-Getty & Veron’s (2001) Gatalog-which 
is still very incomplete at high redshifts- contains more than 10 ^ already identified blazars, 
in addition to more than 10"^ quasars and other less energetic AGNs. Actually, GLAST 
mission enhanced sensitivity is expected to uncover more than 10"^ 7 -ray emitting AGNs 
(Gehrels & Michelson 1999). If the actual number of 7 -ray emitting AGNs below the EGRET 
detection threshold is, say, 10'^ (which can be a conservative assumption: approximately 1 
7 -ray blazar per 10 "^ galaxies), they might produce high-latitude sources lacking a clear low- 
frequency counterpart. Even when considering reduced probabilities for microlensing with 
large magnifications (i.e. ~ t/A^, where r is the local optical depth and A the magnification). 
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we still found that a handful of the sources at high latitudes already observed by EGRET, and 
maybe a few hundreds of sources to be observed by GLAST, can be the result of microlensing. 
For instance, considering a random lensing probability of 5 10“^, and a local value of r of 
0{1) (even larger values are possible because the presence of superluminal components in 
blazars -see below-) and amplifications of ~ 100, we could produce by lensing ~ 5 high- 
latitude sources (as we see below) without counterparts at lower energies. Of course, if a 
large number of distant blazars are, say, one order of magnitude below GLAST sensitivity, in 
such a way that they do not need huge amplihcations in order to be detected, the number of 
possible microlensing cases could very much increase. Assuming an amplihcation of a factor 
10, the number of cases could increase to several hundreds in the forthcoming large GLAST 
sample. 


3. Microlensing model for extended sonrces 


We consider a distant and weak y-ray emitting blazar whose GeV flux is well below 
EGRET sensitivity threshold and whose radio flux is at the mJy level, also below the sensi¬ 
tivity of typical all-sky surveys. Assuming a galaxy interposed in the line of sight, the lens 
equation, in the lens plane, is (Ghang & Refsdal 1979, Schneider et ah 1992) 

r - ro - - diag(K -F y, k - 7 )r - cuq = 0, (1) 


where the coordinate system is centered on the lens, with the orientation of the orthonormal 
basis {ci, 62 } chosen to diagonalize the quadrupole matrix; the source is at rg and the image 
position is r. The fourth and hfth terms in Eq. (1) arise from the deflection imprinted 
by the host galaxy as a whole: k and y are the focusing and the shear of the galaxy at 
the point mass position, respectively. The only effect of Uq is to change the unperturbed 
source position rg by a constant. We shall ignore cuq assuming a source position s = tq -|- cuq 
in the lens plane. Defining new appropriate coordinates in the lens and the source plane 
(Schneider et al. 1992), X = [a/|1 — k + y|/i?£;]r, and Y = [1/R£;a/[W-'k”+^]s, the 
lens equation becomes Y = e ( diag(A, 1) ) X — X/ |X|^, where e = sign (1 — k -|- y) and 
A = (1 — K — y)/(l — K -f y). The solution of the latter equation can be found reducing the 
problem to a fourth order equation for (Schneider et al. 1992). We hnally obtain the 
magnihcation, A = lohs/h, as 


AiXf + Xl)^ + eil-A) (Xf-Xi)-l 


( 2 ) 



5 


For extended circular sources, the magnification is given by (e.g. Han et al. 2000): 

^ _ /o"" /o"' ^)Mr, ip)rdrdip 

where (r, (p) are polar coordinates in a reference frame placed at the center of the source, Vg 
is the radius of the source, X(r, p) is the surface intensity distribution of the source, and Aq 
is the sum of the magnifications of all images. We define R = t/Re and suppose that the 
lens is moving with constant velocity v.^ 

The radius of the internal regions of AGNs depends on the energy as a power law, 
r.y oc (Blandford and Levinson 1995). We then define a reference source with radius 
r^ef and y-ray energy such that R^{E) = i?ref(-E/i?ref)^, where capital letters stand for 
normalized quantities using the Einstein radius. We assume that the intensity of the source 
is uniform, and that its spectrum follows approximately a power law Io{E) = Ij.ef{E/E^ei)~^ 
with ^ G (1.7, 2.7), where Rei is the intensity of the reference source (Krolik 1999). The 
surface intensity distribution of the source will then be Tq{E) = Iq{E)/t[R^{EY, from where 
the magnification can be computed. 

Using that A = ///q, we define J = ///ref = A{E/, the intensity in units of Aef. 
We adopt p = 1.1 , .^ = 2, and a reference source with dimensionless radius /?ref = 2 x 10“^ 
and y-ray energy /7ref = 0.1 GeV (Blandford and Levinson 1995). In Figure 2 we show just an 
example of the results of our numerical computations. The y-ray emission of the background 
AGN can be amplified significantly. This magnification can make an otherwise unnoticed 
source to exceed the detection threshold. Note also that the y-sphere corresponding to 10 
GeV, whose size is similar to the optical emitting regions of some typical AGNs, is negligibly 
amplified, while the lower energy curves -well within the EGRET range- all show magnifi¬ 
cations in excess of a factor 100. This phenomenon also has a particular spectral signature, 
produced by the differential amplification of the different y-regions. The spectral evolution 
(chromaticity) effect is a change in the spectral slope at medium energies. This break, pre¬ 
dicted only for microlensing events, and its peculiar time evolution -it shifts towards high 
energies as times goes by since transit- can be used as a specific test to differentiate this 
from other phenomena (see Figure 2, right panels). 


"‘We choose the origin of time (t = 0) as the instant of closest approach between the lens and the source. 
Then, if the center of the source is placed at b = (^i, ^2) when t = 0 (projected on the plane of the lens), the 
position of any point of the source with polar coordinates (r, y) is si(t) = bi — vtcosB + r cosy and S 2 {t) = 
62 — vt sin 9 + r sin y, where n = |v| , 0 is the angle between v and ei, 0 < r < and 0 < y < 27r. In units of 
the Einstein radius, 17 = [i?i—Tcos 0 -I-i?cosy]/ R \1 — k -|- 7I, and I2 = [B 2 — T sin6*-|-i?siny]/ a /|1 — k -I- y], 
where T = vI/Re and B 12 = bi^ 2 /RE- When 7 = 0 (no shear), we can take 9 = 0 and B = (0,i3o) without 
loosing generality. 
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4. Time scales and number of events 


The time scale of lensing variability is to = Re/v. For instance: for M = O.IMq, Hq=75 
km s“^ Mpc“^, ;2s = 7, zi = 0.1, f2o=0.2, and n=5000 km s“^, to ~ 170 days, while the half 
width of the peaks shown in Figure 2 is about 9 days. This may result in a variable source 
for separated EGRET viewing periods (typically of ~ 15 days each). If, instead. Mi ~ 5Mq 
and V ~ 1000 km s“^, the half width of the peaks is ~1000 days, and the y-ray source would 
most likely be seen as a steady, non-variable detection. ® 


Following Subramanian & Gopal-Krishna (1991), if the lenses have the same mass, and 
the observing period is At, the expected number of events, in the case of a background source 
moving with velocity v, will be 


A = r 1 + 


2At 

TTUtn 


(4) 


where r is the optical depth and u = h/vE = BjR e is usually taken as 1. The second 
term between the brackets corresponds to the increment of the probability of detecting a 
microlensing event due to the movement of the source (see Romero et ah 1995 and Surpi 
et al. 1996). The optical depth can be defined as the ratio of the surface mass density of 
microlensing matter to the critical mass density Scrit = c^Dqs/^t^GDoxDis. For simplicity, we 
shall assume that the lenses are located at a distance from the center of the galaxy much 
smaller than the core radius. Then, k ~ Sc/Scrit, with the central surface density. The 
high surface mass density associated with the core of normal galaxies along with the usual 
assumption that most of this mass is in the form of compact objects naturally leads to high 
optical depths for microlensing. For instance, in the case of the lensed quasar Q2237-031, 
where four images are well-resolved, lensing models indicate values of r ~ 0.5 (Schneider et 
al 1988), which are corroborated by the detections of microlensing-based optical variability 
with relatively high duty cycles (e.g. Gorrigan et al. 1991). Other lensed sources display 
even higher duty cycles (e.g. Koopmans & de Bruyn 2000). 


For a mass distribution of lenses given by N{M) oc M““ for Mmin < M < 
(Salpeter 1955), the number of expected microlensing events by stars with masses in the range 
(Ml, M 2 ), included in the total mass range (Mmin, Mmax), during At days of observations is 


Ml-M2 


—V D{a,Mi,M2) At 

I.lbMKT -;--- 

c B{a, Mmin, Mmax) 30days 


(5) 


^The velocities chosen are typical for extragalactic objects. For galaxies at redshifts 2 ~ 0.1 or bigger, 
their recession velocities very much exceed the one we are using (by more than order of magnitude). The 
latter are relative velocities, involving small traversal components of the recession speed and proper motions 
of the star within the galaxy and the galaxy itself (with respect to its cluster). A worked example (for BL 
LAC 0846+51) can be found in Nottale (1986), for which he quotes 300 < v < 3000 km/s. 
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where the functions B and D are, respectively, 


and 


B{a, Mmin, 


/ \ 2—a 

/ \ 2—a 

f A^max \ 

( \ 

\ Mq J 

\ Mq J 


2 — a 
for a 7 ^ 2 , 


= In 




Mrr 
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D{a,Mi,M2) 



In 


3 — 2q; 
for a ^ 



for a 




3 

2 ’ 


( 6 ) 


(7) 


and we have assumed a case with ~ 0.9 and 2:1 ~ 0.1 to fix the numerical coefficient. 
Clearly, the total number of events will strongly depend (apart from the expected influence 
of r and k.) on the velocity of the source and the index a, usually taken in the range (2, 3). 

A particular interesting case of source candidates are blazars where the high-energy 
emission is produced in a superluminal component with apparent velocity n > c in the 
lens plane (Gopal-Krishna & Subramanian 1991). They will produce 7 -ray sources with the 
highest levels of variability. We hnd that during the EGRET’s lifetime, hundreds of events 
can be expected for optical depths in the range 0.2-0.4. Even for values of r as low as 10“^, 
if the source is apparently superluminal, then the number of expected events in the EGRET 
observing time is above 10. Instead, sources whose velocities in the lens plane are much 
smaller than c (say, about 5000 km s“^) will produce only a few events. A complete study 
of the distribution of lens masses within this model will be presented elsewhere. 


5. Concluding remarks 

In summary, we have found that gravitational microlensing of the innermost regions 
of distant AGNs can produce unidentihed sources compatible with those observed at high 
galactic latitudes. While large amplihcation factors are obtained for 7 -rays, a negligible 
magnihcation results in the radio band. In the case of the optical emission, if it is the same 
particle population giving rise to gamma-rays and optical emission then we can also expect 



magnification of the optical luminosity. The gamma-ray spectral evolution provides a spe- 
cihc signature for the microlensing events that can be used to differentiate this from other 
kind of phenomena. Higher (lensing) optical depth, or the presence of shear, will lead to a 
diversity of light curves. We remark that we have used the Chang and Refsdal (1979) model 
as the gravitational lensing scenario, but that galaxies with denser cores would require a 
more detailed treatment. This would necessarily include the study of caustic patterns of all 
stars at a time in order to get the light curves. We shall explore the lensing model in much 
more detail in a subsequent publication (Eiroa et ah 2002, in preparation). 
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Variability index I 


Fig. 1.— Left; Distribution of the photon spectral index of all unidentihed sources with 
galactic latitudes |6| > 10 deg. A Gaussian £t would give a mean equal to 2.6 and deviation 
equal to 0.3. Most of the sources have steep spectra, consistent with the AGN population. 
Right: Histogram for the variability index I of each of the HL-UnidS (/ = 1 is the mean value 
of I for pulsars). A Gaussian curve would have in this case a mean at /=2, a value already 
2a above the expected one for a non-variable population, and a deviation equal to 0.7. More 
than half of these sources are classified as likely variable. Taking only those sources having 
|6| > 40°, both mean values are even higher. See Torres et al. (2001a,b) for details. 
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Fig. 2 .— Lensing results for a dimensionless impact parameter equal to B/Rg = b/vg = 0.5 
(upper panel) and 2 (lower panel). Notice the specific feature (a changing break in the 
power-law spectrum) predicted for the spectral evolution. In the left panels we show the 
light curves, from top to bottom E= 0.1 GeV, 0.5 GeV, 1 GeV, 10 GeV. In the right ones 
the spectrum evolution, from top to bottom T = vt/RE =0, 0.02, 0.05, and 0.10. Lensing 
parameters were chosen as k = 0.4, 7 = 0 . 







